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Granular solid hydrodynamics, constructed to describe quasi-elastic and plastic motion of granular 
solid, is shown also capable of accounting for the rheology of granular dense flow. This makes it a 
unified, though still qualitative, hydrodynamic description, enabling one to tackle fluidization and 
jamming, the hysteretic transition between elasto-plastic motion and uniform dense flow. 
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A potentially catastrophic feature of granular media is 
their variable capability to sustain external stresses. As 
the mechanical stability of any structure relies on this ca- 
pability, it is important to have a thorough understand- 
ing when and why it weakens, and how it recovers. The 
transition from a solid-like response to a liquid-like one is 
fluidization; and jamming denotes increasingly often the 
reverse change. After a mud slide starts, after it fiuidizes, 
jamming is when it stops again - at what stress and den- 
sity, and whether the village downhill will be reached, are 
then questions of considerable interests. 

The two limiting states of the transition may be re- 
ferred to as granular solid and uniform dense flow. In 
the first, the grains are deformed and at rest, with all 
energy being elastic. In the second, they jiggle, rattle, 
move macroscopic distances, and a significant portion of 
the energy is kinetic. The transition from the dominance 
of one energy to the other may be gradual or abrupt, and 
has two possible intermediate states: the uniform, ideally 
plastic critical state 0, Q , or shear band, the nonuniform 
path. A unified theory of these limiting and intermedi- 
ate states does not as yet exist, though GSH (for granu- 
lar solid hydrodynamics) seems close. It was originally 
constructed to account for granular solid and its elasto- 
plastic motion. Here, we demonstrate its applicability 
to dense flow: GSH displays broad agreement with three 
existing theories on various aspects of dense flow, and 
accounts for the data of Savage and Sayed |3[ . 

The first of the three existing theories is by Pouliquen 
et al. [i| • Starting from the insight that granular rheology 
in dense flows is controlled by a dimensionless parameter 
~ j/y/P (where 7 is the shear rate, P the pressure), they 
distilled two locally applicable constitutive relations from 
experiments and simulations, for the density p and the 
friction angle a s jP (with o s the shear stress), 
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where, with p cp the closed-packed density, p r = p/p cp is 
the relative one. pi,p2 denote, respectively, the friction 
angle for 7 — > 0, 00. The authors took n — 1, though the 



difference to n = 2 is subtle, as both describe a gentle 
change from pi to P2 with 7. 

Earlier, Boquet et al. Q developed a continuum theory 
to account for their experiment. Starting from the results 
of the kinetic theory for inelastic hard spheres, they mod- 
ified the density dependence of the pressure P, viscosity 
77 and relaxation rate 7 to accommodate the higher den- 
sity in their system. They employ the Cauchy stress cr^ , 
and a balance equation for the granular temperature Tg' 



P 



P- 

T G 



Wij, 



dT G /dt ~ 



1 



Pr 



7 Tg , with 
~ (l-fr)"' 



(2) 
(3) 



where a — 1, (3 between 1 and 2.5. «y = i(V \vj + VjVi) 
denotes the strain rate, with its traceless part, and 
v 2 = v%v% the scalarized shear rate. (The notation 
and s are used also for other tensors below, such as the 
strain and stress.) Solving both equations for a Couette 
cell, the solution was found to agree well with their data. 
This theory II does not consider any elastic contributions. 

Considering shallow flows on an inclined plane and ro- 
tating drums, Aranson and Tsimring identified the hys- 
teresis of transition, or the delay between jamming and 
fluidization, as a key feature of granular behavior [f|. 
Their theory III treats the Cauchy stress as the sum of 
two parts, a solid-like, possibly elastic contribution ga^, 
and a rate-dependent fluid one. A crucial variable is 
an order parameter g that is 1 for granular solid, and 
for dense flow. The authors take the friction angle 
(j), differently than above, as the ratio of the solid stress 
components, and postulate a free energy f{g) such that 
granular solid, g = 1, is unstable for large shear stresses, 
<t> > 0i ; while dense flow, g = 0, is unstable for small 
ones, 4> < 4>o- Both are stable in the intermediate region, 
4>\ > tfi > tpo . This theory does not consider variations in 
the density p, or in Tg, and takes afj as an input from 
some other theory. But its success provides a pivotal 
insight: The viability, even appropriateness, of using a 
partially bistable energy to account for the hysteresis. 

GSH starts from the basic fact that grains with endur- 
ing contacts are elastically deformed. Its essential idea 
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is that this deformation is slowly lost when grains jiggle, 
as they briefly loose contact with one another. Granular 
solid's complex elasto-plastic behavior was shown to be 
a result of this simple physics, assuming the dominance 
of the elastic energy. Kinetic energy, or granular heat, is 
what underlies the behavior of granular gas. So it seems 
obvious that the behavior of dense flow results when both 
energies are comparable, when the contribution to the 
stress from granular temperature becomes equally im- 
portant as that from deformation. 

GSH was first employed to calculate static stress distri- 
bution for various geometries, including sand piles, silos, 
and point load, achieving results in agreement with ob- 
servation Q. It was then employed to consider slowly 
strained granular solid, and found to yield response en- 
velopes similar to those from modern hypoplastic the- 
ory [8|. Recently, the critical state - generally consid- 
ered a hallmark of granular behavior - was identified as 
a steady-state, elastic solution of GSH (Qj : Although given 
as a simple analytic expression, the solution realistically 
renders the critical state and the approach to it, including 
dilatancy and contractancy. Finally, the velocity of elas- 
tic waves were calculated as a function of the stress , 
and found to agree well with experiments 

GSH consists of • conservation laws for the energy w, 
mass p, and momentum pvi, • an evolution equation for 
the elastic strain «y, and • balance equations for two 
entropy densities, s and s g . Two entropies are neces- 
sary, because granular media display a two-stage irre- 
versibility: Macroscopic energy, kinetic and elastic, dissi- 
pates into mesoscopic, inter-granular degrees of freedom, 
mainly granular jiggling and the collision-induced, fluctu- 
ating elastic deformation. After a characteristic time, the 
energy degrades further into microscopic, inner-granular 
degrees of freedom, especially phonons. The granular 
and the true entropy, s g ,s, account respectively for the 
energy of the meso- and microscopic degrees of freedom. 
The elastic strain mj is the portion of the total strain £y 
that deforms the grains and leads to reversible storage 
of elastic energy. The rest-frame energy density wq is a 
function of s g , s, p, Uij (though we shall neglect s, as we 
are not interested in effects such as thermal expansion 
at present). The conjugate variables are: Granular tem- 
perature T g = dwo/dsg, chemical potential p = dwo/dp, 
elastic stress, Try = —dwo/duij, and the gaseous pres- 
sure Pt = p 2 d(wo/p)/dp\ s , . The elastic stress 7Ty 
derives from granular deformation, while Pt is gener- 
ated by granular temperature - similar to the tempera- 
ture generated pressure in a gas. All conjugate variables: 
(Tg, p, 7r.y ;, Pt) are given once wq is. 

In GSH, the Cauchy stress <Tjj [given by momentum 
conservation, d(pvi)/dt + Vj(o~ij + pviVj) — 0] and the 
balance equation for s g are given as 



Tij = (1 - a)TTij + P T Sij 

dsg/dt = (r] g v 2 
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Although a result of general principles, the expression for 
o~ij is, remarkably, a simple sum of the elastic stress, the 
gaseous pressure, and the viscous stress, with r\ g the shear 
viscosity. (Compressional flow is usually negligible. If 
not, one needs to include the bulk viscosity.) For elasto- 
plastic motion, only (1 — a)TTij is important; granular 
gas is well accounted for by Pt% — Vg v iji dense flow 
needs all three terms, a ~ 0.8 is a softening coefficient 
that remains constant for all shear rates considered in 
the present context. (It becomes smaller only for ultra 
low shear rates, in ratcheting or elastic waves). Note 
the similarity of Eq (TJJ to the above cited theories, with 
the difference that theory II ignores 7Ty, and theory III 
takes it as given. In Eq 7 is the relaxation rate of 
s g , accounting for the inelastic collisions that occur when 
grains jiggle. The positive term rj g v 2 = r]g v ij v ij describes 
how grains, being sheared past one another, start to jiggle 
in the process, leading to an increase of s g . From a more 
general point of view, this term describes how the kinetic 
energy dissipates into granular heat. In the stationary 
limit, for ds g /dt = 0, we have T g = v s \Jr)gj~i- (Only a 
uniform T g is considered - more terms exist otherwise.) 
See for derivation and detailed explanation. 

Eqs (|4I5[) hold for the given set of variables, indepen- 
dent of the granular material, or the specific form of Wq. 
Material-specific properties are encoded in u>q(p, u^, s g ), 
also the transport coefficients: r\ g ,~l- We obtain them 
from qualitative consideration, also comparison to ex- 
periments and existing theories. (More puristically, one 
would of course like to obtain them from simulation or 
microscopic calculations.) For dry sand and glass beads, 
a simple energy expression, the sum of the elastic energy 
wi(uij,p) and granular heat W2{s g , p): wq = +u>2, has 
turned out to be quite adequate as a first approximation. 
Then 7Tjj = —dw-\ / dug , Pt ~ p 2 d(w2/ ' p) / 'dp (if one as- 
sumes <C I, see [12J), which is why stress and energy 
contributions are simply linked: If w± dominates, only 
Ttij is important; while Pt hinges on a sufficiently large 
w 2 . We have w x = B(p)y/A [A 2 +U 2 ./£], with A = -u u , 
u 2 = vf-jiqj. For a granular system at rest, w\ is the 
only energy, and the elastic stress is the total stress. The 
mentioned calculation of static stress distributions was 
carried out using w\ [l3| . Granular heat W2 is the lowest 
order expansion in s g , 
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The linear term vanishes because granular jiggling dissi- 
pates and decreases toward zero, implying W2(s g ) is mini- 
mal for s g — 0. Expanding also the transport coefficients, 



Vg = Va+ ViTg, 7 = 7o + JiT g , T g = u s vWti. (7) 

we take rjo <C 771X0, 70 <C JiT g , as is appropriate for 
any T„ typical of elasto-plastic motion and dense flow, 
see Then the last of Eqs holds, for ds g /dt = 0. 



3 



In theory II, Tq is the energy per degree of freedom, 
W2 ~ Tq, while w 2 ~ s 2 ~ T g . We therefore iden- 
tify T g ~ a/Tg, and note the perfect agreement between 
Eqs P15|) and Eqs (I4I5I6I7|) . This is important, because 
the Independence (in contrast to the density dependence 
discussed below) is rather fixed - it is the result of the 
kinetic theory on one hand, and the general considera- 
tion rendered above on the other. (Note taking W2 ~ T g 
would disregard the fact that T g dissipates, and W2 is 
minimal for s g , T g = in an adiabatic system. It is more 
appropriate for ideal gas than the granular one.) 

To consider the density dependence, we focus on 1 — p r 
= \—pj Pep-, which represents a stronger dependency than 
p if the sand is dense, p r w 1. We take 
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with a = i, /3 = § to fit the experimental results of [H 
with respect to the polystyrene beads. The expression 
for Pt derives from b = b {l — p r ) a , cf Eq ©. Taking 
b ~ ln(l — p r ) would yield Pr ~ 1/(1 — p r ) exactly, as 
in Eqs (|3]), but also leads to a divergent s g , see [12j |. As- 
suming a w 0.1 approximates the result, yet avoids the 
problem. The coefficients bo,h\,gi are (material depen- 
dent) numbers. Combining Eq (|4|6I7|8|) and denoting the 
shear rate as 7 = V x v y (hence 7 = \f2v s in simple-shear 
geometry) , we arrive at the final expressions for the pres- 
sure P = <7«/3 and the shear stress a s = (rr^rr^A 1 / 2 . 
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where C x = \abp cp p 2 h 1 /g 1 , C = \p cp ^/h\l2g[. P C ,II C 
denote the rate-independent, elastic contributions, with 
ILc/Pc independent of 1 — p r , see the explanation below. 
The first of Eq © may be written as 1 — p r ~ 
= i/s/Pr- It is the same as Eq (fTJ) if P c 
an understandable mismatch, because the 



is neglected 

consideration of theory I involves inertia and confining 
pressure, but neglects elasticity. The friction angle cr s /P 
is given by II C /P C for 7 — >■ 0, and CJ C\ for 7 — > 00, which 
we may respectively identify as p\,p2- The number n of 
Eq (TT]) is 2 in GSH. As mentioned, the difference to n = 1 
is subtle for the friction angle - but less so if one look 
at the pressure or shear stress individually. That both 
grow wiihjj 2 is in fact a behavior that already Bagnold 



observed 14(. Finally, a note on volume versus pres- 
sure control: Yielding P, <r s for given 7, 1 — p r , Eqs © 
are directly appropriate for experiments performed under 
constant volume. If P is fixed, one uses the first calculate 
p, and rewrite the second as a s — U c = (P — P C )C/C\, 
with a coefficient C/C\ that does not depend on 1 — p r 
- though still on p, a weaker function of P and 7. 

Returning to the exponent of 7, we note a + j3 — 2 
if (<r s — II C )/(P — P c ) is to be independent of 1 — p r ; and 



(3 — a = 1, if Pt ~ (1 — p r )~ 2 ■ Together, they imply a = 
= |, as given above, see Figffl However, for glass 
beads of the same experiment Q, Pt ~ 7 2 /(l — Pr), or 
/3 = a = 1 is more appropriate. In addition, the friction 
angle decreases for increasing 7 here, implying II C /P C > 
C/Ci, without contradicting any general principle. 

Next we discuss the elastic contributions: P c , II C . Ap- 
plying a constant shear rate v s to an elastic body, the 
shear stress will monotonically increase - until the point 
of breakage. Sand is different and can maintain a con- 
stant stress, afy This is the famous critical state 0, Q 
that has, for given density, a unique, rate-independent 
stress value. Employing GSH, this is easy to understand: 
Because elastic deformation Uij is slowly lost if the grains 
jiggle, and because grains indeed jiggle when forced to 
shear past one another, a shear rate v s not only in- 
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FIG. 1: Comparison of GSH to the polystyrene data of Savage 
and Sayed [3|, with Pressure P, shear stress a s , and a s /P 
given as functions of the shear rate 7. The first two figures 
show the 7 2 -dependence, and the third the convergence onto 
the weakly density-dependent, high-rate limit p2 = C/Ci. 
Diamonds, squares and circles are the experimental points at 
the specified densities p. (Data for the largest p are not used, 
because the authors believe they may be plagued by "finite- 
particle-size effects.") The curves render Eqs ([9]), with hi = 
3.1 • I0" 4 v / ab^, 91 = 12f.7 v / a6o 3 , a = 0.1, p cp = 0.64p bli;fc . 
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creases Uij , as in any elastic medium, but also decreases 
it. The critical state is the steady state in which both 
processes balance each other, such that the elastic defor- 
mation remains constant over time, in spite of a finite 
v s . As shown in [9], the stationary solution depends 
on the density, but not on v 8 . The associated elastic 
TTn (ut e , p) = TTj p Sij /3 - 7Tj u& / v s , character- 



stress = ni j K u, kl , 
ized by two scalars, 7r^ 
v s . Its contributions in Eq 



and 7rJ, is also independent of 



are: P c = |(1 - a)-K c u , 
II C = (1 — a)n c s . Although both P c , H c depend on 1 — p r , 
the ratio H c /P c = ^s/ n u does n °t- 

We did not find any independent data on the critical 
state of polystyrene beads, though that from [§| indi- 
cate P c s» 50 Pa, (Tc/Pc ~ 0.25, implying that the softer 
polystyrene beads have a B « 10 5 Pa, while the other co- 
efficients retain their orders of magnitude as given in Q . 
(Note: 7Ty ~ B, and B s» 5 x 10 9 Pa for sand.) 

At lower shear rates, say for v s < 10 s _1 , the rate- 
dependent terms of tr^ are quadratically small, Pt ~ 

T 2 

s 



rjiTgV^ ~ w 2 , and may be neglected. This is the 



reason the total stress is given by the rate-independent 
critical state, cry = (1 — c^tt?-, for a fairly broad range of 
shear rates, and why soil mechanic textbooks emphasize 
the rate-independence of granular behavior. 

We note that fluidization, as considered above, is uni- 
form and continuous, without anything resembling "fail- 
ure" or "yield." Starting from a state of isotropic stress, a 
sheared granular system will approach the critical state, 
in the continuous way as calculated in [9] . The end state 
is, more generally, given by Eqs ((9]), though the differ- 
ence to the critical state is evident only at higher shear 
rates. There is an alternative path that goes through 
an energetic instability, eg. the Coulomb yield contained 
in wiiitij), see jl2| . which sets in when the ratio tt s /ttu 
becomes too large. This transition is discontinuous, non- 
uniform, and shear bands necessarily appear. We shall 
consider it in a forthcoming paper. 

Jamming, the reverse transition - a drop of the shear 
rate v s from a finite value to zero at given stress - is 
necessarily discontinuous. In contrast to the authors of 
theory III, however, we do not believe this instability is 
marked by a lower bound of tt s /itu, as elastic solutions 
are perfectly stable at isotropic stresses, ir s = 0. Rather, 
jamming seems an instability that sets in when the den- 
sity is too high to enable a shear flow v s . Although v s 
is not a state variable, T g is, and we have T g ~ v s [cf. 
Eq 0] for any processes slow enough for Eq ((5|) to have 
reached its stationary limit. The appropriate instability 
must therefore be in W2{s g , p) = s 2 g /2pb. If we substitute 
b for b, we have Pt instead of Pr, 
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long as p r is away from 1, and all results above remain 
valid. For p equal to 



Pjam = Pcp(l ~ 26i/o), 



(11) 



however, the convexity of u>2 with respect to p is lost, 



and no finite value of T g ' 
from the condition: d 2 W2/dp 2 
from d 2 f 2 /dp 2 \ Tg ~ £(p 2 £b)\ Tg = 0, where f 2 



is stable. Pj am is obtained 
= 0, or equivalently, 



T s 



w 2 



-pbT 2 /2. Eq (fTTj) is the result to lowest order 



in bi and a. Note p < Pj am imples a lower bound for 
v s if the pressure is given instead of the density, as a 
smaller v s will imply a larger p, see the first of Eq (JSJ. 
On a plane inclined by the angle <p, the friction angle is 
t&nip = <y s /P, with the angle of repose given by (p r = 
^{pjam)- Since the two terms ~ j 2 are negligible by 
then, the angle of repose is given by the critical angle at 

Pjam 

: tanyv = U c (p jam )/P c (p jam ). This is consistent 
with observation, because the critical angle is necessarily 
smaller than the angle at which Coulomb yield sets in. 
All these statements are independent of the specific form 
of 6, which may possibly prove inappropriate - though 
the case for an instability in b(p) seems watertight. 

Summary: Because GSH is capable of accounting for 
elasto-plastic motion, including the critical state, and 
also for dense flow, fluidization and jamming, we believe 
that this hydrodynamic theory, conventionally based on 
conservation laws and thermodynamics, is a viable can- 
didate for a unified theory of granular media. 



With bi small, we may neglect the correction term as 
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